Synthesis of thin-long heavy nuclei in ternary collisions 
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We illustrate the formation of a thin-long structure of heavy nuclei by three-nucleus simultaneous 
collisions within time-dependent density functional theory The impact parameter dependence for 
such formation is systematically demonstrated through clarifications of the difference between binary 
and ternary collision events. A new method for producing thin-long heavy nuclei in the laboratory 
is suggested, as well as the possible formation of the thin- long structure in hot dense matter such 
as that encountered in core collapse supernovae. 
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Generally, synthesis of superheavy elements and many 
other exotic nuclei comes from collisions between two sta- 
ble nuclei, whereas simultaneous collisions among more 
than two stable nuclei, which can occur in principle, have 
not been taken seriously. In experiments using accelera- 
tors, it seems difficult to construct a setup to make three 
nuclei collide simultaneously. Nevertheless it is reason- 
able to expect that a fixed target can be bombarded by 
two beams moving in the opposite directions. In astro- 
physical circumstances, on the other hand, there seem 
many chances to see three-body reactions such as triple- 
alpha reactions The triple-alpha reactions, if occur- 
ring through two-body resonances as in the Hoyle pic- 
ture, are similar to radioactive ion beam experiments [1| 
in the sense that two reactions occur non-simultaneously 
in terms of a time scale of the strong interaction. Even in 
stars and supernova cores, however, simultaneous three- 
nucleus collisions within the strong interaction time scale 
are usually presumed to be rather hard to encounter, let 
alone more than three nucleus collisions. How hard they 
are remains to be examined thoroughly, although the pos- 
sible influence of simultaneous triple-alpha fusion reac- 
tions on carbon production in stars has attracted much 
attention [3[. 

Two-nucleus fusion reactions in stars and supernova 
cores are often investigated in two steps 0, Ha]. One 
starts with fundamental two-nucleus fusion reactions in 
vacuum and then incorporates medium (plasma) effects. 
For example, in the famous Gamow approach, one solves 
the elementary tunneling problem for given relative ki- 
netic energy and then obtain the temperature-dependent 
reaction rate by averaging the resulting kinetic-energy- 
dependent reaction rate over the Boltzmann distribu- 
tion of the kinetic energy. In this work, we concentrate 
on the fundamental binary and ternary collisions within 
the time-dependent density functional theory (TDDFT). 
Among multi- nucleus collisions, ternary collisions are 
particularly worth investigating in the sense that there 
always exists one plane containing all the three center- 
of-mass coordinates of colliding pairs. Such a geometric 



restriction facilitates the formation of a low-dimensional 
quantum system. As we shall see, a thin-long structure 
of the fusion products can be stabilized by rotation in 
the case of non-central binary and ternary collisions. 

The TDDFT approach, which was originally proposed 
by P. A. M. Dirac [6j, is useful for describing nuclear 
collisions at low energies except the reactions below 
the Coulomb barrier. In the present case, we perform 
three-dimensional TDDFT calculations with Skyrme- 
type effective nucleon-nucleon interactions (SLy6 [7[ and 
SKI3 [8]). One can then derive various quantities from 
the many-nucleon wave function self-consistently ob- 
tained in the form of the Slater determinant. We re- 
mark that the colliding system would break into a few 
fragments with various neutron/proton ratios [si for the 
collision energy per nucleon of order the nucleon Fermi 
energy, and multi-fragmentation takes place for rather 
higher energies. Multi-fragmentation, in which two- 
nucleon collisions are usually expected to play an im- 
portant role [10], is out of our scope. 

Since 56 Fe can be abundant in dense stellar matter, we 
consider heavy element synthesis due to multiple 56 Fe col- 
lisions within ~ 33.3 x 10 -22 s, which corresponds to the 
typical duration time of low-energy heavy- ion reactions. 
We first show an example of ternary collisions with the 
initial condition shown in Fig. [T] and the incident energy 
of 1.47 MeV per nucleon in the center-of-mass frame, 
which is illustrated in Fig. [2] In this case, the fusion 
product is 168 Pt. In contrast with reactions producing 
light elements, the corresponding reactions are endother- 
mic and hence require additional energy even after the 
contact. For the impact parameters that allow Nuclei 
I and II as well as Nuclei II and III to touch with each 
other, corresponding to the cases of \b\ = 1, 3, 5, 7 fm, the 
final products have a thin-long structure of length 25-30 
fm, which is stabilized by rotation. For example, in the 
case of \b\ = 7 fm, rotation up to tt rad is achieved in 
a few 10 -21 s, which corresponds to the total rotational 
energy of order 2.10 MeV. Note that thin- long nuclei can 
occur at optimal values of the incident energy, below and 
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FIG. 1: (Color online) The initial positions of three identi- 
cal 56 Fe nuclei, Nucleus I, Nucleus II, and Nucleus III, which 
are set to (|6|,0, 15 fm), (0,0,0), and (-|6|, 0,-15 fm), re- 
spectively. The velocity vector of each nucleus is given as 
(0,0, |v|), (0,0,0), and (0,0, -|v|), respectively. The volume 
of a box inside which the simulation is performed is 48 x 48 x 48 
fm 3 . 

above which fusion is drastically suppressed. Note also 
that the stabilization of such a thin-long structure by 
rotation can be seen in exotic clustering of light nuclei 

0. 

It has been shown here that thin-long structures are 
produced by non-central ternary collisions for a broad 
range of the impact parameter and that the rotation 
speed becomes faster for larger \b\ as it should. Note 
that thin-long structures turn into more spherical shapes 
within 33.3 x 10~ 22 s if there are no rotations, which 
suggests that the rotation plays an indispensable role in 
keeping the fusion products thin-long. 

For comparison, we consider two-nucleus collisions by 
starting with the same initial configurations as in Fig. [T] 
except for the absence of Nucleus II. The corresponding 
incident energy is 2.20 MeV per nucleon in the center- 
of-mass frame, and the fusion product is 112 Te. We can 
observe from Fig. [3] that a similar but shorter structure 
occurs in the case of two- nucleus collisions. The degree of 
deformation from a spherical case is significantly larger 
in the case of ternary collisions (see Table H]). We have 
thus quantitatively confirmed that ternary collisions are 
more efficient than binary collisions in producing thin- 
long structures. 

Note that fusion is not only one of the most efficient 
ways of producing heavy elements in stars, but also the 
main method for producing superheavy nuclei in the lab- 
oratory. The present TDDFT calculations suggest that 
fusion reactions of two or three identical 56 Fe, which lead 
to the synthesis of chemical elements heavier than iron, 
can be easier to occur with the help of rotations. 

The efficiency of multi-nucleus collisions in the produc- 
tion of thin-long heavy nuclei can be essentially under- 
stood by the competition between the surface tension and 
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(b) \b\ = 3 fm 
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(d) \b\ = 7 fm 
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FIG. 2: (Color online) Ternary collisions with different im- 
pact parameters that start with the initial condition shown 
in Fig. [T] (the SLy6 parameter set is taken). The velocity 
is set to be 1.95 x 10 22 fm/s. Snapshots at 1.33xl0~ 22 s 
(left), 13.3 xlO -22 s (middle), and 26.6 xlO -22 s (right) are 
shown for \b\ — 1, 3, 5, 7 fm, while those at 1.33xl0 -22 s, 
13.3xl0~ 22 s, and 20.0xl0~ 22 s are shown for \b\ = 9 fm. 

the Coulomb repulsion. Let us describe this competition 
within the framework of an incompressible liquid-drop 
model, which can characterize a specific geometry of the 
early stage of multi- nucleus collisions (see Fig. HI (a)). In 
this model we restrict ourselves to the situation in which 
the total volume of the system is fixed to ^irr^A/S with 
ro = 1.20 fm and the number of nucleons A. 

Let us now consider a string of n identical nuclei of 
radius R = ro^/n) 1 / 3 , which we regard as a precursor 
of a thin- long heavy nucleus as long as n > 1 (Fig. HI (a)). 
Then, we obtain the length R z of this string as 

R z = 2nR = 2r n 2/3 A 1/s . (1) 

The liquid-drop model allows one to quantify the com- 
petition between the surface tension and the Coulomb 
repulsion, which is essential to the formation of thin-long 
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(a) \b\ = 1 fm 




(b) |b| = 3 fm 




(c) \b\ = 5 fm 




FIG. 3: (Color online) Two- nucleus collisions with different 
impact parameters. The difference from Fig. [2] lies solely 
in the absence of Nucleus II in the initial configurations 



illustrated by Fig. [T] Snapshots at 3.33x10 



(left), 



6.66 xlO -22 s (middle), 13.3 xlO -22 s (right) are shown for 
161 = 1,3,5 fm. 



heavy nuclei. The surface and Coulomb energies of the 
string can be written as 



^Surf = a. 
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A 5 / 3 , (3) 



is 



where k = A/Z with the number of protons Z, 
the distance between the centers of the z-th and j-th nu- 
clei that constitute the string, and a s and a c are set to 
17.0 MeV and 1.38 MeV fm, respectively, in a manner 
that is consistent with empirical masses of stable nuclei. 



TABLE I: Length of the thin-long products by ternary colli- 
sion at 26.6xl0 -22 s and by binary collision at 13.3xl0 -22 s 
that are shown in Figs. [2] and [3] The diameter for a min- 
imum sphere that contains each product is calculated and 
then averaged over all the cases in which the product is fused. 
The resultant length, denoted by 7Z, is divided by A 1 ^ 3 . For 
comparison, a typical diameter for a spherical nucleus is also 
given. 
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FIG. 4: (Color online) An incompressible liquid-drop model 
for multi-nucleus collisions, (a) Linearly connected n spher- 
ical nuclei with the total and proton numbers (A,Z). (b) 
d n E(K, = 2, n, A) is depicted as a function of n (A = 162, 199, 
212, 225, 262). A state with n = 1 loses its global stability 
at A = 225, which is a little bit larger than Ao(k = 2) = 212, 
and, for A > 225, a state with n > 2 is stabilized. 



Accordingly, the competition is characterized by the sum 

E(K, n, A) EE £ S urf + ^Coul 



17.0 n^ 3 A 2 ^ + 



2.30 



-2/3 



+ (2- n - l)n 



-5/3 



A 5 / 3 (MeV). (4) 



Although n and A are integers, for optimization we take 
n and A as real numbers. The extremal condition for 
given k and A can then be obtained from 

d n E(K,n,A) 

850^ 2 n 2 2 n - [(184n - 575)2 n + 239n + 575] A 



150ft 2 n 8 / 3 2™A- 2 / 3 



(MeV). 



The mass number satisfying d n E(K, n, A cr i t ) 
responds to 



t (ft,n) 



850^ 2 n 2 2 n 



(184n - 575)2™ + 239n + 575 ' 



(5) 
cor- 

(6) 



TABLE II: The energy E(k, n, A) calculated from Eq. (gj) as 
function of n. We take A — 400, 500 and k — 2.5, which are 
typical values expected for ternary collisions of heavy nuclei 
in the laboratory. 





n=3 


71=2 


71=1 


£(2.5,n,400) 
£7(2.5, n, 500) 


3.28 GeV 
4.38 GeV 


3.30 GeV 
4.45 GeV 


3.32 GeV 
4.55 GeV 



where the denominator is positive for any n. At least one 
set of the numbers (ft, n,A) satisfying d n E(K,n, A) = 
exists only when A > A CY i t (K,n). Otherwise, 
d n E(K,n, A) > holds. Since, for given ft, A cr i t (ft, n) 
has the global minimum value Ao(ft) = 2 12 (ft/2) 2 , we 
can conclude that the value of Ao(ft) gives a rough cri- 
terion for the occurrence of the stability transition from 
the state with n = 1 to the state with n > 2. In fact, for 
A < Aq(k) where d n E(tv, n, A) > is satisfied for any n, 
a single sphere (n = 1) is energetically preferred, while 
for A > Ao(ft), as shown in Fig. HJb), there can be more 
than one value of n that fulfills d n E(K,n, A) = 0. 

As a result of detailed analyses of the energy landscape, 
we find that a state with n = 1 remains energetically 
optimal for mass number up to a value that is a little bit 
larger than Ao(ft). Above this value, the optimal state is 
not preferably realized in two-nucleus collisions (n = 2) 
but in multi- nucleus collisions (n > 3). This implies that 
for superheavy synthesis, production of thin-long nuclei 
by ternary collisions (n = 3) can be more efficient than 
the production by usual binary collisions leading to the 
states with n = 1 or n = 2 (Table HI)). 

In summary we have found a thin-long structure of 
heavy nuclei as a result of simultaneous three-nucleus 
collisions within the TDDFT approach. The validity of 
the calculations using the SLy6 interaction has been con- 
firmed by using the other interaction (SKI3 [8]) in terms 
of the realization and rotational stabilization of thin-long 
heavy nuclei. This study is expected to provide a motiva- 
tion for designing a new accelerator and detector system 
for superheavy synthesis in which three-nucleus simulta- 
neous collisions can take place. In fact, in addition to 
the existing methods for superheavy synthesis that are 
based on "binary" fusion and multi-nucleon transfer re- 
actions, "ternary" fusion reactions and subsequent rota- 
tional stabilization of compound nuclei could provide a 
novel method for the superheavy science. 

It is also interesting to consider possible astrophysical 
implications. For example, the three- nucleus fusion rate 
in stars and supernova cores can be estimated by allow- 
ing for the average over the initial configurations through 
the Boltzmann factor and the plasma effects, i.e., electron 
screening effects and many-body Coulomb correlation ef- 
fects between ionic nuclei Q. The latter effects act to 
reduce the Coulomb barrier between the colliding nuclei 
and thus to enhance the fusion rate. If the plasma is rela- 
tively dilute and hot as in the Sun, the plasma effects can 
be safely ignored. In this case, an extension of the usual 
Gamow rate to three- nucleus fusion applies, leading to 
the fusion rate as a function of the plasma temperature 
T. Once the plasma becomes dense and thus strongly 
coupled as in supernova cores [H|, the plasma effects 
should manifest themselves in the fusion rate through 



the three-particle static correlation function for electron- 
screened ions. However, it is a challenging problem to 
accurately obtain the three-particle correlation function 
as a function of T and the plasma density p [IK . Recall 
that there are optimal values of the incident energy for 
three-nucleus fusion in vacuum. In matter, the optimal 
values can be lowered by the plasma effects, leading to 
enhancement of the three-nucleus fusion rate through the 
Boltzmann factor. If the inverse of this rate at the val- 
ues of p and T relevant for stars and supernova cores is 
shorter than the corresponding evolutionary time scale, 
one can expect that thin-long heavy nuclei occur in such 
celestial objects. 
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